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ABSTRACT 

We propose a method to estimate the pseudo jet opening half-angle of GRBs using 
the spectral peak energy (ifpj-peak luminosity relation (so called Yonetoku relation) 
as well as the .E p -collimation-corrected 7 -ray energy relation (so called Ghirlanda re¬ 
lation). For bursts with known jet break times and redshifts, we compared the pseudo 
jet opening half-angle with the standard one and found that the differences are within 
a factor 2. We apply the method to 689 long GRBS. We found that the distribution 
function of the pseudo jet opening half-angle obeys f(9j) oc Qj 2 - 2±q - 2 with possible 
cutoffs for 9j < 0.04 and 9j > 0.3 although the log-normal fit is also possible. 9~ 2 dis¬ 
tribution is compatible with the structured jet model. From the distribution function 
we found that the beaming correction for the rate of GRBs is ~ 340, which means 
~ 10 -5 yr -1 galaxy -1 or only one in 10 2 type Ib/c supernovae. We also found the 
evolution of the distribution function as a function of the redshift. 
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1 INTRODUCTION 

Gamma-ray bursts (GRBs) arise from re lativistic jets (e.g., 
iMeszarosI 120021 : IZhane fc Meszarod l2003lf . Relativistic mo¬ 
tion, with Lorentz factor of greater than ~ 10 2 is necessary 
in order to resolve the “Compactness problem”. On the other 
hand, the evidence of the jet collimation has been derived 
by the wave-length-independent achromatic break in the ob¬ 
served afterglow light curve. The jet opening half-angle has 
been actually measured in this context for bursts with suc¬ 
cessfully measured redshifts and jet break times. However, 
the small number of samples prevents us from discussing 
the statistical properties of the opening angle distribution. 
We have not yet known the maximum and minimum value 
of the opening ang le as we ll as the slope of the distribu¬ 
tion function well llFrail et alJl200D . These quantities are 
very important to argue the GRB rate, the energetics, the 
nature of progenitors, and so on. 

There is a correlation between the rest-frame spec¬ 
tral peak energy E p and the isotropic equival ent 7 -ray en ¬ 
ergy Ej so of GRBs called the Amati relation jAmatLet_a]J 
120021 : lAtteiall2003l : lLamb et, alJl2Qo3 : ISakamotoeUffTl^oiMlh 

Several authors have argued theore t ical interpretations 
of the relation IlYamaz aki et al. ]20M jEjchler fc Levinsonl 
120041 : lLamb et~alTl2005k While NakarA^piranlll20(uii ~ar-^ 
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gued against the Amati relation because ~40% of 751 
BATSE GRBs do not have the solution to the undeter¬ 
mined redshift under the Amati relation and clear outliers 
such as GRB 980425 and GRB 031203 exist. A similar 
ar gument against t he Amati relation has a lso bee n made 
bvlBandW^PreecjjIgOOSlb jGhirlanda et alJ (|200Ml_Jj£SHgd 

agains t l^aka^^Piranl 11200411 as well as iBanrT fc Preec3 

ll2005ll . They used 442 bright BATSE GRBs with good 
stati stics a nd pseudo redshift from the lag-luminosity rela¬ 
tion iBand. Norris fc Bonnell2004lh and obtained the Amati 
relation with the slightly different power law index and the 
larger scatter than the original one. They found that the 
chance probability of the revised Amati rela tion is 2.1 x 
10 -65 . Similar analysis can be also seen in Bosniak et al. 

B . 

As for outliers, lYamazaki. Yonetoku fc Nakamm5 
i2003f) argued that if the intrinsic E p of GRB 980425 is 
2-4 MeV similar to GRB990123 and the viewing angle 
is ~ 9j + IO 7 -1 with 9j and 7 being the jet opening 
half-angle and the gamma factor of the jet, the observed 
E iso and E p of GRB 980425 can be reproduced. If we locate 
GRB 980425 at the distance larger than 800 Mpc, the event 
is not observed. Hence, GRB 980425-like events are not 
included in the bright BATSE GRBs. Sim ilar arguments 
have been also done for GRB 031203 llR.amirez-R.uiz et all 
l2004ll . Note here that in the classification of HETE-2, 
GRB 980425 and GRB 031203 are X-ray rich GRBs. 
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Recently, similar relations to t he Amati one w i th the 
tighter correlation have been found. lYonetoku et all (2004) 
have shown that for 16 GRBs with known redshifts detected 
by BATSE and BeppoSAX, there is a tight positive corre¬ 
lation between E p and the peak luminosity ( E p -L p rela- 
tion). The chance proba bility is extremely low 5.3 x 10~ 9 . 
iGhirlanda et alJ ( 2005btl called the relation as the Yone¬ 
toku relation and checked the validity of the relation using 
442 bright GRBs with the pseudo redshift and confirmed 
the relation with the same power law index within the er¬ 
ror in the original Yonetoku relation. They found that the 
chance proba bili ty of the Yonetoku relation is 1.6 x 10 -69 . 
lYonetoku et alJ (2004}) used the Yonetoku relation as the 
redshift indicator for 689 GRB samples without known red- 
shift, and derived the GRB formation rate. On the other 
hand, for 15 GRBs detected by BATSE, BeppoSAX and 
HET E-2 with m e asure d redshift and opening half-angle 9j , 
iGhirlanda. et alJ ( 2004 1 found that E p correlates with the 
collimation-corrected 7 -ray energy E-y. We call the relation 
as the Ghirlanda relation. 

The Ep-Lp diagram of GRBs may correspond to the 
Hertzsprung-Russell (HR) diagram of stars. Then the Yo¬ 
netoku as well as the Ghirlanda relations correspond to the 
main-sequence where stars cluster around a single curve de¬ 
termined by the mass of the star. This suggests the existence 
of a certain parameter that controls GRBs like the mass of 
the star in the HR diagram. In the HR diagram the out¬ 
liers of the main sequence exist such as red giants and white 
dwarfs. We know the physical reasons for the existence of 
these outliers so that they do not refute the main sequence 
relation. In the Yonetoku and the Ghirlanda relations also 
outliers exist for which we will know the phy sical reasons 
in future as suggested in (lYama, zaki. Yonetoku fc Nakamural 
120031 iRamirez-Ruiz et alJl2004th 

In this letter, we show that the Yonetoku and the 
Ghirlanda relations can be used to estimate the opening 
half-angle of the relativistic jet of GRBs. Main advantage 
of our jet opening half-angle estimator is that the open¬ 
ing angle can be calculated only from the information of 
the prompt emission. This letter is organized as follows. I 11 
section [ 5 ] using the Yonetoku and Ghirlanda relations, we 
derive an empirical formula to estimate the opening half¬ 
angle and discuss its validity. Using the estimator, we derive 
the distribution of opening half-angle of BATSE-triggered 
bursts in section [3 Section 0] is devoted to discussions. 


2 JET OPENING ANGLE INFERRED FROM 
TH YONETOKU AND THE GHIRLANDA 
RELATIONS 

Usually the jet opening half-angle is estimated only when 
both the redshift and the achromatic break time in the after¬ 
glow light curve are measured. Under the simple assumption 
of a uniform ambient matter distribution of number density 
no, the jet opening half-angle is estimated as 

break = 0.12 [tjet,d/(l + 2)] 3 ^ 8 (^0»?7/E'i S o,53) 1 ^ 8 , (1) 

where Ei BO ,53 = Ei BO / 10 53 erg, and tj e t,d and r ; 7 are the 
jet break time in days and the efficiency of the fireball in 
converting the en ergy in the ejecta into 7 -rays, respectively 
llSari et all Il99it) . However continuous follow-up observa¬ 


tions are required to measure the achromatic jet-break time, 
and moreover these kinds of observations are realized only 
for the bright afterglows so that it is hard to measure the 
opening half-angle for large amount of GRBs in this method. 

We propose here a different method to estimate the 
opening half-angle using only informations of the prompt 
gamma-ray emission. Let us assume that the rest-frame 
spectral peak energy E p , the peak luminosity L p , the 
isotropic equivalent 7 -ray energy E \ so , the jet opening-half 
angle 9j , and the collimation-corrected energy E 7 satisfy the 
Yonetoku and the Ghirlanda relations as 

E p = 2.1 x 10 2 L p , 52 °' 50 keV , (2) 

E p = 4.8 x 10 2 E 7 , 6 i °- 71 keV , (3) 

E 1 = (1 — cos9j)E iso , (4) 

wh ere Q x d enot es Q/ 10 x in cgs units (lYonetoku et al.ll2004l : 
IGhirlanda et aljl2004 ). From equations (2), (3) and (4), we 
have, 

= 1 - cos 9j = 0.30 Lp, 52°- 50/0 ' 71 £- o 1 i51 • (5) 

From the observed E p and L p , we can determine the redshift 
under the Yonetoku relation and the given cosmological pa¬ 
rameter as in (lYonetoku et ali2004l) . Then Ei so can be com¬ 
puted using the observed fluence. Therefore, only from the 
information of the prompt gamma-ray emission, we can esti¬ 
mate the jet opening half-angle for each GRB. This method 
has a strong advantage compared with the jet-break mea¬ 
surement in afterglows since we can use the large number of 
BATSE GRBs data. 

To show the validity of our method, in figure 0 we 
compare 9j, break estimated by Eq. © with 9j by Eq. 0 
for GRBs with measured redshift and jet break time. It is 
found that there exists a positive correlation with the lin¬ 
ear correlation coefficient including weighting factor is 0.760 
with 13 degree of freedom (excluding the lower and upper 
limit samples), which corresponds to the chance probability 
of 9.63 x 10 -4 . Their differences from the linear function 
(equivalent line) are within a factor of 2. Although we had 
better use the word “the pseudo jet opening half-angle” for 
the half-angle from equation (4), we simply use the opening 
half-angle in this letter for convenience. 


3 JET OPENING ANGLE DISTRIBUTION OF 
BATSE BURSTS 

We used the same data set of 689 long GRBs published in 
( Yonetoku et alJ|2004T ) to measure the distribution function, 
f(9j). Here 9j is the jet opening half-angle and f{9j)d9j is 
the relative number of the jet with the opening half-angle 
between 9j and 9j + d9j . In order to have a better signal- 
to-noise ratio in our analyses, we selected 605 GRBs with 
the flux greater than 1 x 10 -6 erg cm - 2 s -1 . Having obtained 
the observed peak flux and the spectral indices as well as the 
E p , the redshif t is estim ated using the Yonetoku relation for 
each event llYonetoku et a UEiill:. Then, we can estimate 
L p and Ei so . Hence 9j is also calculated by equation 0. 
The empirical 9j estimator shown in equation 0 depends 
on two parameters; L p and Ei so . This fact means that the 
two different types of selection effects are mixed. To avoid 
such doubly truncation effects, we manually set a truncation 
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Figure 1. Comparison between jet opening half-angle Oj esti¬ 
mated by Eq. 10 and 0j, break estimated by Eq. 0. The solid 
line is an equivalent line. 



Figure 2. Distribution of the jet opening half-angle Oj estimated 
by Eq. 0 vs. redshift derived from the E p -L p relation. The solid 
line shows the truncation of the upper bound of Oj that is caused 
by the flux limit of Fu m = 2 x 10 —7 ergs cm“ 2 s _1 . The cross 
points are the Oj measured in the present work. The opening 
angles measured from the jet-break time in the optical afterglow 
are also plot on the same figure as squares and arrows. 

as the Eiso > L p x 1 sec at the rest frame of GRBs. In this 
case, both L p and Ei so are limited only by the peak flux, so 
we can deal with the data as the simply truncated one. 

Figure 0 shows the distribution of the derived jet open¬ 
ing half-angle as a function of the redshift. The Oj dis¬ 
tributes within the range of 0.04 < Oj < 0.3 radian. Using 
the flux limited samples explained above, we can estimate 
the opening half-angle evolution and the true distribution 
of Oj following the detailed mathem a tical descriptions (e.g. 
Lvnden-Belll Il97ll: IPetrosianl 1 19931: iMalonev &; PetrosiarJ 
^99Tf^netokiT etalll2004f) . 

First, we estimate the opening half-angle evolution, 
which is the redshift dependence of f(0j ), from the (1 + 
z,9j) distribution shown in figure [5] The data correla¬ 
tion degree in the flux limited samples are calculated by 


Figure 3. The best index of the jet opening angle evolution 
measured by the r-statistical method. When we assume Oj evo¬ 
lution as the function of (1 + z) k , k = —0.45 is the best value. 
We also show 3 sigma upper and lower bound of the index as 
k = —0.26 and k = —0.63, respectively. No evolution of 0j was 
rejected about 6 a confidence. 


so called r-statistical method which is very similar to 
the Kendell’s r statistics. To refer to the previo us works 
iLlovd-Ronning et alJ I20 qT lYonetoku et, al.lf2004l . we as¬ 
sume the functional form of f(0j) evolution as (1 + z ) k , 
and calculate the data correlation degree for each k value. 
In figure [3 we show the (1 + z,6j) correlation as a func¬ 
tion of index k. No 6j evolution is rejected about 6 sigma 
confidence level. When we assume the case of k = —0.45, 
Oj becomes independent of the redshift. In other words, 0j 
evolution of (1 + z)~ 0A5 is hidden in the (1 + z, 0j) plane. 

Next, we derive the opening half-angle distribution. 
On e possib le method (non-parametric method) developed 
bv lLvnden-Belll lll97ll) is applied in this analysis. We define 
Oj = 0j/{ 1 + a )” 0,45 which is equivalent to the evolution- 
removed opening half-angle. Then the observed data ran¬ 
domly distributes in the (1 + z,0j) plane, so we can easily 
assume the missing data caused by the flux limit. Strictly 
speaking, for high-redshift GRBs (e.g., a > 5), the observed 
fluence may be estimated as the lower value because there 
may be missing photons behind the background level. On 
the other hand, the peak flux may be correctly estimated. 
Hence derived opening half-angles for high-z events may be 
larger than the actual values. In order to avoid such confu¬ 
sion, we consider 430 events within the estimated redshifts 
of less than 4.5. Additional merit is that the k-correction fac¬ 
tor can be neglected. The differential Oj distribution /(#'•) 
which corresponds to one at z = 1 , is shown in figure 21 If 
one would like to obtain the 0j distribution at each redshift, 
it is roughly estimated as f(0j) = /(#')(( 1 + z)/ 2)~ 0 ' 45 . 

The opening half-angle distribution /(#'■) is based on 
the number of event detected by BATSE. The chance prob¬ 
ability to observe narrowly collimated events are much lower 
than that of the wide opening half-angle events due to the 
geometrical effect. Therefore, we have to take into account 
the correction of f(0'j)/£lj when we argue the true probabil¬ 
ity density function. In figure 0 we show the true opening 
half-angle distribution. It can be fitted by the power-law 
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ing the same method with smaller number of samples. Since 
the lower cut-off of the distribution 9'j ~ 0.04 are sim¬ 
ilar for both cases, the diffe rence comes from the higher 
power law index ~ —4.54 of iFrail et all (l200lh . Since the 
number of GRBs is much smaller in their analysis, obser¬ 
vational biases of the smaller samples of larger opening 
half-angle causes the steeper power law^ Adopting the lo¬ 
cal GRB rate of ~ 0.5 Gpc _3 yr _1 iSchmiddl200lll . we ob¬ 
tain the true rate of GRB as ~ 170 Gpc _3 yr _1 which means 
~ 10~ 5 yr _1 galaxy _1 with the number density of the galaxy 
being ~ 10 -2 galaxy Mpc -3 . Since the rate of type Ib/c su - 
pernovae is ~ 10~ 3 yr _1 galaxy _1 dGappellaro et a I] 12003b 
only one in 10 2 type Ib/c supernovae becomes GRB if GRB 
is the peculiar type Ib/c supernova. 


Figure 4. The differential 0'. distribution at 2 = 1 of the GRBs 
observed with BATSE detector. 



Opening Angle [radian] 

Figure 5. True opening angle distribution at z = 1. When we 
adopt the power-law model for the region 0j > 0.05 radian, we 
obtain the best result of g'y 2 - 2 ± 0.2 w j|p gg ty 0 confidence level. 

form in the range 9'j > 0.05 radian as g'- 2 - 2±0 - 2 , Clear cut 
off at ~ 0.Q4 rad can be seen . 

iGhirlanda et alj (f2005a l obtained the opening half¬ 
angle distribution using the Ghirlanda relation and the Am- 
ati relation. The result looks like a log-normal one. We also 
tried to fit a log-normal one. We obtained acceptable results 
with the mean log/x = —1.15 (0.07 radian) and the stan¬ 
dard deviation logo- = 0.24 (0.04-0.12 radian). However the 
data distribution is asymmetric in the logarithmic horizontal 
scale as shown in Figure^] To say the power law distribution 
with cut-offs is better than the log-normal one, we have to 
confirm that the cut-off at 9j = 0.04 exists. However we can 
not conclude the existence because of the large errors in the 
small 9j region. 

From the distribution function we found the beam¬ 
ing correction for the rate of GRBs ~ 340 at a = 1 
by computing = f f(9j)Q^ 1 d9j / f f(9'j)d9j. This 

rate is about factor of 4 larger than that obtained by 
lGuetta 1 Piran fe Waxman 1 2004 ) using a different method, 
while obtained this correction ~ 520 us- 


4 DISCUSSIONS 

We present an empirical opening half-angle estimator that is 
inferred from the Yonetoku and the Ghirlanda relations. Our 
method requires only the data of prompt emission, which is 
different from the standard method using the break in the 
afterglow light curve and the redshift. Using the empirical 
opening angle estimator, we have derived the opening angle 
distribution of GRBs. The distribution can be fitted by the 
power-law form in the range 9j > 0.07 rad as gj 2 - 2±0 - 2 . The 
cut off at ~ 0.04 rad can be seen. In the uniform jet model 
of GRBs, this means that the distribution function happens 
to be ~ 9~ 2 . 

The other possibility is the structured jet model. In 
the original version of the universal structured jet m odel 
llRossi. Lazzati fc R.eeslliool : IZhang fc Meszarojl2003l . the 

energy per unit solid angle is in proportion to ~ 9~ 2 . The 
lower and upper angle cutoff exist and the jet structure is 
essentially the same between the lower and the upper cut¬ 
off, which is the origin of the name “structured” jet. The 
viewing angle corresponds to the jet opening angle in the 
structured jet model so that we need to argue observation- 
ally what will be the opening a ngle distribution of the struc¬ 
tured jet iPerna. Sari fc Fraill 1200.^1 . If all bursts were ob¬ 
servable, the distribution would be uniform per unit solid 
angle and f(9) oc 9. However, Ei BO for the smaller viewing 
angle is brighter by a factor of £U 2 so that the maximum 
observable distance is larger by a factor of 0 _1 which con¬ 
tains a volume larger than 9~ 3 . Then we have f(9) oc 9~ 2 
which is compatible with the result in th e present paper. 
Espec ially our Fig. 5 looks like Fig. 3 of lPerna. Sari fc Fraill 
(2003) with appropriate parameters. 

The evolution effect found in the present paper means 
that the jet opening half-angle becomes narrower for larger 
redshift. One possible qualitative explanation is that this is 
due to the metal dependence of the progenitors. Since the 
metalicity of the star decreases as a function of the redshift, 
it is expected that the mass loss of the stars decreases as 
a function of the redshift if the mass loss is derived by the 
line absorption of photons in the atmosphere of the stars. 
This suggests that for high redshift progenitors of GRBs, 
the mass of the envelope is larger so that only the sharper 
jet can punch a hole in the envelope of the progenitor star. 

In this paper, when we calculate the true opening half¬ 
angle distribution, we implicitly assumed that the jet emis¬ 
sion can be seen only when the jet is seen on-axis. In re- 
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ality, effects of ofF-axis viewing G RBs might be important 
(lYamazaki et alJl2002L 1200 3. l2004jh . When the contributions 
of off-axis emission is considered, true distribution may be 
modified. However, the beaming correction from off-axis ef¬ 
fects is importa nt only for z < 1.5 bec ause of the relativistic 
beaming effect liYa.ma.za.ki et a.lJl2004li . 
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